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Chapter  1 


Introduction 

Narrowband  and  wideband  ambiguity  functions  are  fundamental  in  the 
study  of  parameter  estimation  by  active  radar  and  sonar  systems.  Interest  in 
ambiguity  function  theory  has  been  renewed  due  to  the  current  intense  research 
activity  in  afSne  wavelet  transforms  and  other  group  theoretic  transforms.  We 
have  attempted  to  continue  the  program  of  unifying  scattered  results  in  radar  and 
sonar  ambiguity  theory  via  group  represent' tion  theory.  The  underlying  idea  is 
that  both  the  narrowband  and  wideband  ambiguity  functions  are  coefficients  of  the 
unitary  representations  of  their  respective  groups.  Wideband  ambiguity  functions 
are  coefficients  of  the  affine  group  and  narrowband  ambiguity  functions  are  coeffi¬ 
cients  of  the  Heisenberg  group.  This  fact  provides  insight  into  important  concepts 
of  admissible  signals,  ambiguity  conservation,  and  ambiguity  function  invariance 
properties  that  are  important  for  signal  design. 

The  history  of  the  applications  of  group  theory  to  radar  and  sonar  is  not 
long  since  the  subject  has  been  less  than  fifty  years  old.  The  stage  opened  in  1953 
with  P.  M.  Woodward’s  seminal  book.  Probability  and  Information  Theory,  with 
Applications  to  Radar[AA,  pages  115-125],  in  which  he  detailed  the  radar  ambiguity 
function  and  its  uses.  Volume  conservation  was  noted,  and  the  chapter  closed  with 
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the  problem  of  signal  design,  based  on  the  desired  radar  ambiguity  function  shape; 
the  radar  synthesis  problem.  Wilcox  [43]  attacked  the  problem  via  the  techniques  of 
mathematical  analysis  with  some  success  in  1960,  but  it  was  not  until  around  1983 
that  Walter  Schempp  started  publishing  a  series  of  papers  culminating  in  his  article 
[39j  and  book  [38].  These  papers  detailed  the  connection  of  the  radar  ambiguity 
function  with  the  Heisenberg  group  and  attacked  the  problem  from  the  viewpoint 
that  the  radar  autoambiguity  function  is  a  positive  definite  function  on  its  group. 
A  lot  was  known  about  positive  definite  functions  already,  especially  through  the 
work  of  Naimark  and  Gelfand,  summarized  in  the  former  author’s  work  [31].  In 
1985,  L.  Auslander  and  R.  Tolimieri  [4]  took  Wilcox’s  aforementioned  paper  and 
studied  it  from  the  group  theoretical  viewpoint  [4]  (as  suggested  by  Schempp,  since 
they  have  cited  one  of  his  preprints?).  Since  the  Heisenberg  group  had  a  discrete 
subgroup,  they  could  apply  nilmanifold  theory,  a  theory  they  had  written  about 
extensively  earlier  [3]. 

About  1985,  Sibul  and  Titlebaum  started  investigating  the  wideband  ambi¬ 
guity  function  [37],  following  the  pioneering  work  by  R.  Altes.  In  1987,  Sibul  with 
the  help  of  R.  Urwin  wrote  a  proposal  for  funding  in  investigating  the  role  of  group 
theory  with  signal  processing.  After  finding  that  the  affine  group  was  the  fun¬ 
damental  concept  for  the  wideband  ambiguity  function,  we  found  that  Auslander 
and  Gertner  detailed  the  same  idea  and  more  in  [2].  Furthermore,  Auslander  and 
Gertner  described  the  approximation  of  the  wideband  to  the  narrowband  in  such  a 
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clear  way  that  we  began  to  wonder  if  there  was  any  connection  between  the  affine 
and  the  Heisenberg  groups.  By  this  time,  A.  Banyaga  joined  us  as  mathematical 
advisor.  In  1989,  we  decided  that  one  had  to  be  a  contraction  of  another.  But 
within  minutes  after  our  meeting,  Banyaga  found  a  preprint  of  Kalnins  and  Miller 
[26]  in  the  lounge  showing  that  our  suspicion  was  indeed  right.  However,  we  did 
find  some  major  conceptual  errors,  although  the  basic  idea  is  correct.  Their  article 
has  been  published,  still  with  those  errors.  So  part  of  this  thesis  is  to  correct  them. 

As  this  thesis  was  being  written,  there  were  ongoing  questions  in  the  signal 
processing  community.  One  of  those  questions  was  concerning  the  existence  of  at 
least  two  forms  of  wideband  ambiguity  functions  and  when  and  how  to  use  them 
properly.  We  found  that  these  ambiguity  functions  arose  from  different  forms  of 
the  affine  group,  and  so  have  included  the  result. 

Finally  we  include  some  results  of  a  another  investigation  in  the  applications 
of  Group  theory.  We  wanted  to  follow  the  methodology'  of  Auslander  and  Tolimieri 
[3]  in  their  treatment  of  the  FFT  to  create  a  transform  based  on  the  finite  affine 
group. 

A  brief  overview  of  the  thesis  now  follows:  Chapter  2  will  sketch  some  ideas 
concerning  the  narrowband  and  wideband  ambiguity  function.  Chapter  3  then 
reviews  some  mathematics  necessary  for  the  rest  of  the  thesis  and  then  shows  how 
the  ambiguity  functions  are  related  to  group  theory.  We  conclude  the  chapter  with 
the  unification  of  the  various  forms  of  the  wideband  ambiguity  in  the  engineering 


4 


literature.  Chapter  4  shows  how  the  wideband  ambiguity  function  may  approximate 
the  narrowband  under  suitable  conditions,  and  then  shows  how  contractions  may 
explain  the  phenomenon.  Chapter  5  then  applies  a  celebrated  theorem  to  ambiguity 
volume  and  then  shows  how  ambiguity  volume  depends  on  contractions.  Finally, 
our  results  on  the  finite  affine  group  are  stated  in  the  appendix. 

Our  original  contributions  consist  of  corrections  of  Kalnins  and  Miller’s  re¬ 
sults  concerning  contractions.  The  second  is  the  unifying  different  forms  of  the 
wideband  ambiguity  function  in  the  engineering  literature.  The  third  is  the  in¬ 
vestigation  of  ambiguity  volume  with  contractions.  And  the  final  is  an  explicit 
construction  of  the  group  theoretic  transform  based  on  the  finite  affine  group. 


Chapter  2 


The  Ambiguity  Functions 


2.1  Introduction 

This  chapter  provides  a  sketch  of  ambiguity  function  theory.  It  is  not 
intended  to  be  exhaustive,  but  to  be  a  motivation  for  the  rest  of  the  thesis.  A  good 
reference  for  further  information  is  (12,  30,  45]. 

2.2  Narrowband  and  Wideband  Echoes 

Radar  and  active  sonar  use  the  same  principle.  The  systems  consist  of  a 
transmitter/receiver  with  a  processing  unit.  A  waveform,  called  a  signal  is  trans¬ 
mitted  towards  an  object  of  interest,  such  as  an  airplane  preparing  to  land  on  the 
runway.  After  a  time  delay,  the  transmitter/receiver  receives  a  reflected  form  of 
the  signal,  called  the  echo,  from  which  the  processing  system  may  extract  impor¬ 
tant  information  such  as  the  position  and  radial  velocity  of  the  object.  Radar  uses 
electromagnetic  waves  in  the  atmosphere,  while  active  sonar  uses  acoustic  waves  in 
water. 

It  is  assumed  that  the  environment  of  the  system  is  free  of  clutter.  The 
waveform  travels  from  the  system  to  the  object  and  back  without  any  interference. 
Furthermore,  we  assume  the  object  to  be  a  point.  The  reflected  waveform  will  be 
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identical  to  the  transmitted  waveform,  delayed  of  course,  if  the  object  were  not 
moving.  Finally  we  assume  the  object  to  be  travelling  at  a  constant  radial  velocity, 
V.  The  ramifications  of  this  assumption  will  be  seen  in  section  4.2. 

Finally,  a  signal  s(f)  may  be  seen  essentially  as  an  element  of  L^(7Z).  It  heis 
time  as  its  domain  and  voltage  as  its  range.  The  signal's  Fourier  transform  may  be 
interpreted  as  the  set  of  its  frequencies,  or  spectrum. 

The  three  assumptions  above  provide  a  basis  for  radar/sonar  theory  in 
general.  Suppose  that  the  signal  was  narrowband  ais  well.  Narrowband  signals  are 
those  whose  Fourier  transforms  have  very  small  support,  called  the  bandwidth{BW), 
so  that  they  essentially  appear  to  be  impulse  functions  in  the  frequency  domain. 
This  also  means  that  their  Fourier  transforms  are  concentrated  about  a  central 
frequency.  We  also  assume  that  the  signal  has  a  very  short  time  duration  T.  The 
narrowband  assumption  is  coupled  with  the  time  duration  assumption  as  follows 
[41,  page  241]: 


£W<_2_ 

c  TBW’ 


(2.1) 


where  c  denotes  the  speed  of  the  signal  in  the  signal  in  the  medium.  Under  this 
additional  assumption,  the  echo  takes  the  form 


e„(<)  =  3{t  +  T)c"*". 


(2.2) 


Here,  r  denotes  the  time  elapsed  between  the  transmission  of  the  signal  and  re- 
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ception  of  the  echo,  called  the  time  delay,  from  which  the  position  relative  to  the 
transmitter/receiver  may  be  calculated.  Let  c  denote  the  speed  of  the  signal  in  the 
medium.  Then  the  position  R  of  the  object  is  calculated  as: 

(2.3) 

since  we  assumed  that  the  object  is  travelling  at  speeds  much  less  than  c.  We  divide 
by  2,  taking  account  that  the  signal  had  to  travel  to  the  object  and  back. 

The  symbol  u;  denotes  the  doppler  shift  and  basically  is  related  to  the  radial 
velocity  of  the  object.  For  example,  the  pitch  of  a  car  moving  towards  an  observer 
will  be  perceived  higher  that  when  it  is  moving  away.  The  radial  velocity  of  the 
object  given  the  doppler  shift  may  be  calculated  using  the  formula: 

2v 

Wo - =  W,  (2.4) 

c 

where  wo  is  the  essential  frequency  of  the  transmitted  signal. 

Thus  we  may  process  the  return  signal  to  find  the  object’s  radial  velocity 
and  range. 

Without  the  additional  narrowband  assumption,  the  return  signal  takes  the 

form: 


c«,(t)  =  \/as{at  +  6) 


(2.5) 
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where  6  is  related  io  the  delay  of  the  first  transmitted  photon  and 


a 


1±£ 

1-0 


(2.6) 


where  0  =  Note  that  now  the  time  scale  of  the  return  signal  is  dilated  or  con¬ 
tracted  depending  on  whether  the  object  is  moving  away  or  towards  our  transmit¬ 
ter/receiver.  In  fact,  nature  is  wideband,  meaning  that  the  narrowband  formulation 
is  an  approximation.  We  will  explore  this  further  in  section  4.2. 


2.2.1  The  Ambiguity  F\inctions 

A  final  assumption  is  that  a  correlation  receiver  produces  forms  a  sufficient 
statistic  for  detection  in  real  white  Gaussian  noise  [40,  page  269].  The  location 
and  velocity  of  the  object  may  be  determined  by  “comparing”  the  echo  with  shifts 
of  the  transmitted  signal.  The  maximum  value  corresponds  to  the  location  and 
velocity  of  the  object. 

The  correlation  between  two  signals  /(<)  and  g{t)  is: 


r  (2.7) 

J^oo 

recognizable  as  the  inner  product  of 

Basically,  to  determine  the  location  and  velocity  of  the  object,  one  simply 
correlates  the  echo  with  shifts  of  the  transmitted  signal: 
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Ce,,(a,T)=/  e{i)s{at  +  T)dt.  (2.8) 

J^^OO 

The  parameters  (a,  t)  that  maximize  the  value  of  Ce.,  are  those  which  best  models 
the  received  signal.  The  location  and  speed  of  the  target  can  be  readily  calculated. 
Each  correlation  process  is  similar  for  each  type  of  echo. 

Abstracting  our  correlation  process,  we  obtain  the  narrowband  and  the 
wideband  ambiguity  functions  Nj^g  and  Wjy. 


Na^.y)  =  r  (2.9) 

=  r  +  b)dt  (2.10) 

oo 

These  two  functions  are  the  key  to  many  synthesis  and  design  problems  in 


radar. 


Chapter  3 


Group  Theory  and  Radar  and  Sonar 

3.1  Introduction 

After  introducing  the  signal  ambiguity  functions  in  the  previous  chapter, 
we  are  now  in  a  position  to  see  how  group  theory  provides  a  unified  view.  We  review 
some  mathematical  background,  with  special  emphasis  on  relevant  techniques,  and 
then  show  how  signal  ambiguity  functions  arise  as  connected  to  special  groups. 
Finally,  we  unify  the  different  forms  of  the  wideband  ambiguity  functions. 

3.2  Some  Mathematical  Background 

Basic  knowledge  of  group  theory  and  differentiable  manifolds,  along  with 
the  theory  of  Lie  groups  and  Lie  algebras  will  be  assumed.  We  also  use  knowledge 
of  elementary  Hilbert  space  freely.  However,  for  clarity,  we  will  review  the  concepts 
of  group  representation  theory  and  also  the  little  known  topic  of  parametrization. 

The  set  of  operators  on  a  Hilbert  space  that  preserve  its  inner  product 
forms  the  unitary  group  under  the  operation  of  composition.  Specifically,  the  space 
will  be  our  Hilbert  space  under  consideration. 

It  is  possible  to  define  measures  on  topological  groups  and  hence  Lie  groups, 
such  as  the  Heisenberg  group  and  the  affine  group.  These  measures  are  compatible 
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with  the  group  structure  in  the  sense  that  if  a  set  is  translated  by  multiplication  of 
a  group  element,  its  measure  remains  unchangeu.  There  can  be  left  invariant  and 
right  invariant  measures,  depending  on  whether  the  set  is  multiplied  by  the  group 
element  on  the  left  or  on  the  right.  Explicitly,  if  /  is  a  measurable  function  on  a 
group  G,  the  measure  /i  is  left- invariant,  if 


figo*9)Mg)  =  /(5)d/i(</),  (3.1 ) 

for  all  go  €  G.  The  right  invariant  measure  is  defined  similarly.  If  an  invariant 
measure,  left  of  right,  exists  for  a  group,  it  is  well  known  that  it  must  necessarily 
be  unique  up  to  multiplication  by  a  constant.  Finally,  if  the  left  and  right  invariant 
measure  coincide,  the  group  in  question  is  said  to  be  unimodular. 

A  unitary  representation  U  of  a.  group  G  on  a  Hilbert  space  H  Is  a  mapping 
assigning  to  each  group  element  a  unitary  operator  U{x)  on  H,  such  that  for  any 
x,y  €  G, 


U{x*y)  =  U[x)oU{y)  (3.2) 

(Recall  that  U  is  a  unitary  operator  on  H  if,  for  all  v,w  €  H,  <  Uv,Uw  >=< 
v,w>  ). 

A  subspace  S  of  H  is  invariant  under  U  if  for  any  u  €  5,  t/(x)  •  u  €  5  for  all 
X  £  G.  A  representation  U  is  said  to  he  irreducible  if  the  only  invariant  subspaces 
of  H  are  the  zero  subspace  and  H  itself. 
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G'  X  G'  — ►  G' 

^  Ip  X  Ip  I  0 

G  X  G  G 

Figure  o.l:  Calculation  of  new  group  operation  under  new  coordinates. 

The  coefficient  of  a  continuous  unitary  group  representation  f/  of  G  on  a 
Hilbert  space  H  with  respect  to  the  ordered  pair  {f,g),  both  /  and  g  in  H,  is  a 
mapping  Cuj,g  :  G  — ►  Tl  : 

X' — *<f,U{x)g>.  (3.3) 

The  coefficient  is  a  generalization  of  the  usual  Fourier  coefficients. 

Finally,  we  cover  the  area  of  coordinates.  Given  a  group,  with  an  operation, 
a  one-to-one  transformation  rp  from  the  group  onto  itself  will  change  its  operation. 
Refering  to  figure  refcomdiagram,  the  group  G  and  G'  are  the  group  with  the 
original  operation  *  and  the  new  group  with  the  new  operation  •  respectively. 
Beginning  with  ★,  we  may  calculate  •  by  simply  noting  that: 

•  =  tpo*o{ipx  ip)~'  (3.4) 

This  technique  is  essential  for  the  correct  contractions  on  the  group  level. 


as  will  be  shown  later.  Corwin  and  Greenleaf  [13,  pages  14-16]  is  one  of  the  few 
references  detailing  this  technique. 
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3.3  The  Affine  and  Heisenberg  Groups 


The  subgroup  of  G/(2,R)  consisting  of  matrices  of  the  form  I  I ,  with 

\0  1/ 

a  >  0  and  b  any  real  number,  is  called  the  real  (two  dimensional)  affine  group. 

It  will  be  denoted  by  A\.  Its  identity  element  is  the  usual  identity  element  of 

/  a  b\  ( a“*  — a~*6\ 

G/(2,R),  and  the  inverse  of  an  element  I  I  I  •  Sometimes 

Vo  1/  V  0  1  ; 

U  b\ 

we  may  denote  an  element  I  I  by  just  (a,  6).  We  mav  thus  see  the  affine 

\o  ly 

group  as  the  set  R+  x  R  with  the  operation; 


(a,b)*(a\b')  =  (aa\ab'-hb). 


We  also  will  need  to  work  with  the  affine  group  cross  the  real  line,  G^,  the 


group  of  matrices  of  the  form: 


a  b  0 


0  10, 


,  0  0  I 


where  again  a  >  0,  6  and  x  are  any  real  number. 

Another  matrix  group  of  particular  interest  is  the  so-called  Heisenberg 
group,  H,  consisting  of  matrices  of  the  form: 


1  X  z 


0  1  y 
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Clearly,  the  identity  and  the  inverse  of  this  set  belongs  to  Ti..  Here,  is  a  three 
dimensional  manifold.  Alternately,  we  may  see  H,  as  with  the  operation 

(i,  y,  z)  *  (x',  y\  z')  =  (x  +  x',  y  +  y',  z  +  z' -h  xy').  (3.8) 

However,  H  has  another  form,  if  one  lets  the  0  of  the  equation  3.4  be 

(x,y,2)i — ►  (x,y,2- ixy),  (3.9) 

we  obtain  a  new  multiplication  rule; 


(x,y,2)*  (x',y',2')  =  (x  +  x',y +  y',2  +  r' +  \{xy'  -  x'y)),  (3.10) 

which  is  called  the  symplectic  parametrization  of  7i.  Needless  to  say,  there  are 
many  other  forms  of  the  Heisenberg  group. 

Both  groups  play  a  fundamental  role  in  signal  processing  via  their  unitary 
representations  on 

Proposition  3.3.1  The  affine  and  Heisenberg  groups  have  unitary  representations 
on  L^'R)  as  follows: 


a.  :>I,-W(m)) 
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U^A<^,b)f(t)  =  ^f(—)  (3.11) 

y/a  a 

b.  VG^-.GA-*U{LHn)) 

UaAaMm  =  ^e''“/(— )  (3.12) 

y/a  a 

c.  U-H:H-^U{C^(n)) 

Un{x.yyZ)f{t)  =  +  x)  (3.13) 

Proof 

The  proof  consists  in  verifying  the  definition  of  group  representations: 

Proof  of  a.: 


(o'. »')/«)  = 

y/a'  a' 


y/^'  \  aa'  ) 


=  £/^.(na',fl6'  +  6)/(0. 


(3.14) 

(3.15) 

(3.16) 


(3.17) 
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Proof  of  b.: 


t/G^(a,6,x)C/G^(a',6',i0/(O  = 


t/o^(<.,6,x)(^e’'“7(!-^))  (3.18) 

(t—b  _  u\ 

—  (3.19) 

V(l^) 


1 


.2irix 

:e  e 


^2iti(x+x' 


=  UG^(aa\ab'  +  b,x  +  2-')f(t). 


(3.21) 


Proof  of  c.: 


i/w(x,y,2)t/wx',y',57(/)  =  C;„(x,y,x)e2-*'c2-V«/(<  +  a.^)  (3.22) 

=  +  x  +  x')  (3.23) 

_  g2iri(*+z'+xv')e2ir,(»+v')«y(^  +  X  +  x')  (3.24) 

=  t/w(x  +  x',y +  y',x  +  c'  +  xy')/(0-  (3.25) 


We  may  now  see  that  the  wideband  and  narrowband  ambiguity  functions 
are  simply  the  coefficients  of  the  affine  and  Heisenberg  groups  respectively. 

But  W/^  does  not  strictly  look  like  an  ambiguity  function.  The  next  section 
will  help  clarify  this. 

Hence  we  may  use  group  theory  as  an  approach  in  analyzing  me  behavior 
of  ambiguity  functions.  In  fact,  Schempp  and  others  studied  the  autoambiguity 
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functions  through  this  view  point,  since  such  functions  are  positive  definite  func¬ 
tions  on  the  groups,  and  many  results  concerning  these  were  available  through  the 
efforts  of  Gelfand  and  Naimark  [31]. 

3.4  Various  Forms  of  the  Wideband  Ambiguity  Function  Unified 

The  form  of  the  wideband  ambiguity  function  above  appears  in  the  guise 
of  wavelets  as  well  in  works  by  Daubechies  and  Miller:  we  shall  call  it  W*: 

^  r  /(()9(— )*.  (3.26) 

yjd  «/— oo  (X 

Speiser  and  others  deal  with  a  second  form  of  the  wideband  ambiguity  function: 

W^{aM  =  P^f{t)g{at  +  b)dt.  (3.27) 

This  is  the  “usual”  wideband  ambiguity  function.  The  third  form: 

W^{x,y)  =  c-5  r  mg{e-t  -  y)dt  (3.28) 

OC 

appears  in  Heil  and  Walnut,  and  subsequently  in  the  thesis  by  Fowler.  And  the 
fourth  form,  which  is  close  to  the  heart  of  engineers  in  signal  processing,  initially 
appeared  in  Kelly-Wishner  and  Altes; 

W*{a,  ^)  =  >/“  /  ^  /(05[a(<  +  ^)1*- 


(3.29) 
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We  would  like  to  show  that  all  these  forms  are  coefficients  of  the  affine 
group,  albeit  with  different  parametrizations.  To  this  purpose,  we  propose  the 
definitions  of  four  groups: 


Definition  3.4.1  Let  the  four  groups  be  defined  as  follows: 

a.  The  matrix  group  Ai  consisting  of  the  set  of  matrices  of  the  form 
with  real  entries  and  s  >  0. 

b.  The  matrix  group  A2  consisting  of  the  set  of  matrices  of  the  form 
with  real  entries  and  s  >  0. 

c.  The  set  with  the  operation: 


S  T 
0  1 
s  0 
T  1 


(x,  y)  *  (x',  y')  s=  (x  +  i',  e~^'y  +  y') 


is  the  group  A3. 

d.  The  set  R+  x  R  with  the  operation 


(3.30) 


(a,6)  *  (a',6')  =  (oa',6+ — )  (3.31) 

a 

is  the  group  Aa- 

Propocition  3.4.1  A1-4  all  isomorphic. 

Proof.  1.  The  group  Ai  is  isomorphic  to  .4^,  by  the  isomorphism  which 


sends  a  matrix  to  its  inverse  transpose.  Indeed,  consider  ^1,2  :  At  -*  A^,  where 
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(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 


The  homomorphism  (t>\,2  is  clearly  one-to-one  and  onto. 

2.  The  group  A2  is  isomorphic  to  ^3  by  the  isomorphism  ^2,3  :  ^2  -^3 

(a  0\ 

where  I  1  *-*  (--lna,6),  and  changing  variables  Ina  •-»  a.  Indeed, 

\b  l) 


^,3 


=  (-ln(aa'),a'fr+6') 

=  (-lna-lna',c'"“'6-|-6') 


=  (— Ina,  6)(— Ina',  6') 


(3.37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 


The  map  02,3  is  certainly  one-to-one  and  onto.  Now,  we  just  change  vari¬ 


ables:  let  X  =  —  Ina  and  x'  =  —  Ina'.  Hence: 


(— Ina,  6)(— lna^60  =  (—Ina  —  lna^a'6^-  6') 

becomes: 


(3.42) 


(x,  6)(x',  b’)  =  (x  +  x',  e-^'b  +  b%  (3.43) 

the  rule  of  A3. 

3.  The  group  A2  is  isomorphic  to  A4  by  the  isomorphism:  (a,  6)  *-♦  (a,a6). 

Indeed, 


02,4(a,6)02.4(a',6')  =  (a,a6)(a',a'6') 

(3.44) 

=  {aa\a'ab+  a’b') 

(3.45) 

=  {aa\a'{ab+ b') 

(3.46) 

=  02.4(aa'ia'(a6+ 6')) 

(3.47) 

Proposition  3.4.2  The  following  maps  are  representations  of  A\  through  A4  on 
L^iH)  respectively: 
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=  i/(— )  (3.48) 

v/a  a 

U2{a,b)f{t)  =  y/af{at-\-b)  (3.49) 

U3(x,y)m  =  e-^fie-‘t-y)  (3.50) 

lUixMit)  =  V^fia[t  +  b])  (3.51) 

Proof  The  proof  consists  of  simple  verification,  Jis  follows: 

a.  The  verification  for  Ui:  done  in  proposition  3.3 

b.  The  verification  that  f/j  is  a  representation  of  A2  is  as  follows: 

{/2(a,6)t/,(a',6')/(0  =  ^^6/2(0, 6)/(a'<  +  6')  (3.52) 

=  y/^'f{a\at  +  6)  +  6')  (3.53) 

=  y/^'f{aa't  +  a'b-\-b')  (3.54) 

=  6/2(00',  a'6  4- 6')/(0  (3.55) 


c.  The  verification  that  U3  is  a  representation  of  A3  is  as  follows: 
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Uz{a,b)U^{a\V)f{t)  =  t/3(a,6)e-t/(e-“'f-6) 

(3.56) 

=  e-^  [e-H  -  b)  -  b) 

(3.57) 

=  c-®^/(e-“+“'<-c-“'6-6') 

(3.58) 

=  f/3(a-Ha'.e-“'6  +  6')/(0 

(3.59) 

d.  The  verification  that  LU  is  a  representation  of  At  is  as  follows: 


U^{a,b)U4{a\y)m  =  C/4(a,6)v^/[a'(t  +  6')]  (3-60) 

=  v^^/[a'{a(/ +  6)  +  6'}]  (3.61) 

=  \/aa'f[a'ai  +  a'ab  +  a'b']  (3.62) 

=  L'  (3.63) 

=  U,{aa',h+-)m  (3.64) 

a 


Hence,  forming  coefficients  as  usual,we  have  all  the  wideband  forms  in  the  engi¬ 
neering  literature  (that  we  know  of). 
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3.5  Invariance  Properties  of  the  Ambiguity  Functions 

We  now  apply  the  concepts  of  the  previous  sections  to  unify  other  results. 
In  [1],  we  find  various  invariance  properties  of  the  wideband  ambiguity  function  as 
detailed  by  Altes.  He  used  the  formula: 

f(t)g[a{t  +  b)]dt,  (3.65) 

OO 

and  so  we  will  reprove  some  of  his  statements  using  the  group  theoretic  viewpoint. 

The  key  idea  is  that  a  unitary  representation  of  a  group  G  preserves  inner 
products: 


<  U{x)f,U{x)g  >=<  f,g  >  .  (3.66) 

The  trick  is  to  be  careful  to  keep  track  of  which  form  one  is  using  as  to  which  group 
one  is  working  with. 

Proposition  3.5.1 

lV;>,6)=lV;,(l/a,-<.6).  (3.67) 

Proof 


=  <f^U,{a,b)g> 


(3.68) 


=  <  U^{a,br^f,U4(a,brW,ia,b)g  >  (3.69) 
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<  i\{lla,-bla)f,g  > 

(3.70) 

<  9,U4{lla,~ab)f  > 

(3.71) 

(3.72) 

What  we  have  done  is  see  complicated  integrals  as  coefficients  of  representations 
of  affine  groups.  We  emphasize  that  the  above  theorem  is  true  for  only  iU  as  the 
inverse  of  an  element  of  another  group  can  take  a  different  form.  For  example, 

W}Ja,  b)  =  Wjj{\la,-bla).  (3.73) 

Some  other  identities  are  special  cases  of  the  following  proposition: 

Proposition  3.5.2 


<  V(x)!,U(y)s  >=<  S,V-'{x)V{y)g  >  .  (3.74) 


Proof  This  proposition  is  a  consequence  of  the  unitary  property  of  the  represen¬ 
tation. 

For  example,  in  the  case  of  the  group  A\ ,  we  obtain  the  results: 


<  U,{a,V)!,V,(ii\b')9  >=<  >, 

a  a 


a 


(3.75) 
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which  may  be  verified  by  the  technique  above.  Thus  group  theory  provides  a  short¬ 
hand  notation  to  bypass  many  integration  substitutions. 

Studying  equation  3.75  further,  we  find  that  some  of  Altes'  identities  arise. 
For  example, 

Proposition  3.5.3  If  f{t)  =  a^u{ai)  and  £f(t)  =  aiv{ai),  then: 

Wl(a,b)  =  WUa,b).  (3.76) 

Proof 

We  note  that  the  hypothesis  implies  that  /  =  f/i(i,0)u  and  that  g  = 
Ui{^,0)v.  Hence: 


<  f,Ui{a,b)g  > 

(3.77) 

<f/,(p  0)u,f/,(a,6)f/,(p0)v> 

(3.78) 

<f/,(p0)«,C/,(p6)u> 

(3.79) 

(3.80) 

<  u,Ui{a,kb)v  > 

(3.81) 

(3.82) 

The  equation  3.80  is  the  result  of  applying  proposition  3.5.2.  Again,  we  must 
emphasize  that  the  form  of  the  equation  will  vary  depending  on  the  group. 


Chapter  4 


Contractions 


4.1  Introduction 

Contractions  form  the  topic  of  this  chapter.  After  motivating  the  idea  of 
contraction,  we  apply  it  to  the  affine  group. 

Contractions  were  presented  for  the  first  time  by  Wigner  and  Inonu  [42]  in 
1953,  to  try  to  connect  special  relativity  with  Galilean  relativity  and  to  connect 
classical  mechanics  with  quantum.  Not  long  after,  Saletan  [35]  provided  a  solid 
foundation  for  Inonu  and  Wigner.  Various  papers  were  then  published  on  contrac¬ 
tions  of  specific  groups,  and  soon  all  the  contractions  of  dimension  three  or  less  were 
known  [11].  In  fact,  the  contraction  we  are  going  to  be  dealing  with  was  listed  in 
that  article  [11].  We  follow  in  particular  the  contraction  definitions  of  Dooley  [15], 
since  they  are  quite  clear  and  sufficiently  rigorous  for  our  purposes.  This  chapter 
contains  corrections  to  Kalnins  and  Miller’s  paper  [26].  Their  errors  lay  primarily 
in  the  improper  parametrization  in  the  Lie  group  level  and  exponentiating  the  Lie 
algebra  under  the  assumption  that  it  was  a  commutative  Lie  algebra,  which  it  is 
not.  After  corrections,  the  contraction  computations  flow  beautifully. 
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4.2  Physical  Motivation 

Sibul  and  Titlebaum  [37],  among  others,  found  that  under  suitable  condi¬ 
tions,  the  wideband  ambiguity  function  may  approximate  the  narrowband.  The 
essential  facts  are  as  follows: 

Let  us  consider  the  wideband  amh’guity  function: 

=  Vs  f  f(t)g{st  +  T)dt.  (4.1) 

Now  let  us  also  assume  that  g{t)  is  a  finite  energy  signal  with  envelope  u{t)  and 
carrier  frequency  wq: 


g{t)  =  u(<)c'^‘  (4.2) 

Hence  the  received  signal  would  have  the  form  of  time  delay  and  dilation: 

g{st  +  t)  :=  u{st  +  T)e'^‘+^>  (4.3) 

Now  recall  that  s  —  where  0  =  ^  .  Expanding  in  Taylor  series: 

s  =  1  -h  2^  +  2/3^  +  •  •  •  (4.4) 

Now  suppose  v  ^  c  so  that  0  4^1  and  ^  1,  then  the  received  signal  becomes 
approximately: 


28 


<7(a<  +  T)  a:  u(<  + 
=  5(<  + 


(4.5) 

(4.6) 


The  associated  ambiguity  function  becomes: 

=  r  /(Oc'“^"p(t  +  r)c‘“^‘dt,  (4.7) 

y—oo 

which  is  the  narrowband  ambiguity  function  up  to  a  phase  factor.  It  is  interesting 
that,  this  approximation  by  physical  considerations  has  a  mathematical  counterpart. 
Essentially,  the  physical  argument  has  shown  that  the  coefficient  of  the  affine  group 
may  approximate  that  of  the  Heisenberg  group.  But  the  key  idea  is  that  the  affine 
group  may  approximate  the  Heisenberg  group  in  a  certain  sense.  This  is  essentially 
the  concept  of  contractions. 

4.3  Contractions 
4.3.1  Definitions 

Definition  4.3.1  Let  Q,Q\  be  two  Lie  algebras  with  the  same  underlying  vector 
space  U.  Then  Qx  is  a  contraction  of  Q,  if  there  exists  a  continuous  mapping 
tfx :  (0, 1]  — »  QC{U)  such  that 
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for  all  x,y^U . 


\im<i>l^[(f>xx,<Psy]c  = 


(4.8) 


Proposition  4.3.1  Suppose  Q  =  {U,  [,))  is  a  Lie  algebra,  with  4>x  as  above  with  A  > 
0.  Then  Qx  =  [»]a)i  where  [i,y]A  =  <f>x^[<i>xi-,4>\y],  is  a  Lie  algebra  isomorphic 

to  g. 


Proof  Since  <t>x  is  invertible. 


(ar,y]A  == 


(4.9) 


is  equivalent  to 


^A[x,y]A  =  [^A3:,</»Ay],  (4.10) 

and  this  certainly  shows  that  ^  is  a  homomorphic  Lie  algebra  to  ^a-  Finally,  <l> 
being  a  member  of  gC{U)  shows  that  it  is  one-to-one  and  onto. 

Proposition  4.3.2  The  map  A  ♦  [,]a  from  (0,1]  to  the  space  of  alternating 
bilinear  forms  on  U,  to  which  the  Lie  brackets  belong  in  particular,  is  a  continuous 
map. 

Proof  This  is  clear  from  seeing  this  map  as  a  composition  of  the  map  <t)x  and  the 
Lie  bracket,  both  of  which  are  continuous. 
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Proposition  4.3.3  Suppose  Q  is  a  Lie  algebra  with  underlying  vector  space  U  and 
(<^a)a€^?+  C  QC{U)  is  such  that 

[i,y]i  =  (4.11) 

exists.  Then  [,]j  is  a  Lie  bracket  on  U  and  the  resulting  Lie  algebra  Q\  is  a  con¬ 
traction  of  Q. 

Proof  We  need  to  show  that  [,]j  is  a  Lie  bracket. 

1.  Linearity. 


[aA'  +  6r,Z],  =  lOT[aA' +  6F,  Z]a 

=  lim  afA',  Z]a  +  limfefK  Z]a 

A— 0  A—O 

=  a^Y,Z],  +  6[r,Z],. 


(4.12) 

(4.13) 

(4.14) 


2.  Antisymmetry 


=  limlx,y'l» 

A— HI 

=  lim-inA'] 


(4.15) 

(4.16) 


=  -(r,Ai,. 


(4.17) 
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3.  Jacobi  Identity 

((A',  y]i,  Z]t  +  [[¥,  Z]u  Xh  +  [[Z,  X],  K],  =  (4.18) 

lim{[[:^,  Y],,  ZU  +  I(A',  YU  Z]x  +  {[X,  YU  Z]x}  =  (4.19) 

0  (4.20) 

Proposition  4.3.4  (Inonu-Wigner  contraction^  LetQ  be  a  Lie  algebra  with  under¬ 
lying  vector  space  U,  and  lei  \  be  a  Lie  subalgebra.  Let  V  denote  a  subspace  of  U 
complementary  to  f.  Thus  x^U  can  be  uniquely  written 

X  =  If  +  xv-  (4.21) 

For  ,  define  a  map  d>\  €  QC{U)  by 

<l>\{x)  =  X|  -f  Axv.  (4.22) 

Let  x,y  be  in  U.  then  (limA—o^r’I^AX,  ^Ay)  exists  and  is  equal  to 


+  [xv',yt]v. 


(4.23) 
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Proof 


(4.24) 

=  lim<^J*[<^A(a:t  +  xv),^A(yt  +  yv)] 

(4.25) 

=  lim^xM®t  +  '^®v,yt  +  <^yv] 

(4.26) 

=  ([®t.  J/tl  +  yt]  +  A[xt,  yv]  +  x^yv,  yv']) 

(4.27) 

But  t  is  a  Lie  si 

ibalgebra,  so  [xt,y|]  €  f  and  so: 

~  (kt^ytl  +  ^((a’l'.ytlt  +  I^fvsytlv) 

(4.28) 

+A(Irtiyv]t  ■+•  (xt,yv'li')  +  A^([xv,yu)t  +  lyv,yv]v)) 

(4.29) 

=  (^t»  ytl  +  l^v',  ytlv  +  yvl^ 

(4.30) 

Note  that  for  the  limit  to  exist,  f  has  to  be  a  Lie  subalgebra.  Otherwise, 

^A*((a^t.yt])  =  ^r(laft^yt]t  +  [3^t»ytlv) 

(4.31) 

=  [a^t^ytlt  +  ^kt.ytl. 

(4.32) 

showing  that  when  X  tends  to  0,  the  limit  does  not  exist. 

Thus  {,]i  is  a  Lie  bracket  on  U,  and  the  resulting  Lie  algebra  is  a  contraction 
of  (?.  A  contraction  arising  this  way  is  called  an  Inonu-Wigntr  contraction  with 
respect  to  f. 
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Definition  4.3.2  Let  G  be  a  Lie  group  and  I\  a  Lie  subgroup,  reductive  in  G.  Then 
Q  may  he  written  as  t0  V,  where  the  decomposition  is  Adfc -invariant  (  Adf^V  C  V 
).  Let  V  M  A'  denote  the  semidirect  product  of  V  by  K  relative  to  this  action.  Define 
Tx  :  V  X  K  G  by  vxivk)  =  expc(Au)it  for  each  A  €  71.'*' . 

Proposition  4.3.5  The  Lie  algebra  ofVtA  K  is  precisely  Q  of  proposition  4-3.3. 
And  the  differential  of  ir,\  at  the  identity  is  precisely  the  map  p\  of  definition  4-3.1. 

Proof  The  first  statement  follows  from  the  fact  that  the  decomposition  is  Adjc- 
invariant.  The  second  follows  from  applying  the  definition  of  differential  of  a  smooth 
map. 

Thus  we  can  define  the  contractions  of  on  a  group  level  as; 

Definition  4.3.3  The  semidirect  product  V  M  A'  is  called  a  contraction  of  G  with 
respect  to  K,  and  the  family  (»’a)a€?i+  of  Tnaps  V  ^  K  G  is  called  the  family  of 
contraction  maps. 

In  practice,  once  we  obtain  a  contraction  on  the  Lie  algebra  level,  we  derive 
a  contraction  of  the  Lie  group  level  simply  by  exponentiating  the  basis  involved. 
We  do  not  know  of  any  general  method  of  contracting  their  coefficients. 

4.4  Contractions  of  the  Affine  Group  to  the  Heisenberg  Group 

4.4.1  Contractions  on  the  Lie  algebra  Level 

We  now  focus  on  Ga  =  xTt,  the  group  of  3  x  3  matrices  of  the  form: 


a  b  0 


0  1  0 


.0  0  e' 


(4.33) 


where  a,b,x  €  '72..  The  group  is  essentially  the  affine  group,  however,  crossed 
with  the  real  line,  which  acts  as  its  center. 

Since  is  a  matrix  group,  its  Lie  algebra  can  be  realized  as  a  matrix 
algebra.  A  set  of  basis  elements  for  Qji  consists  of; 


1  0  O' 


'0  1  O' 


0  0  O' 


=  0  0  0  ,/l2  =  0  0  0  ,>l3  =  0  0  0  ;  (4.34) 


.0  0  0, 


.0  0  0, 


.0  0  1 


with  commutation  relations: 


[Ai,A2]  —  A2,[Ai,A3]  =  [AjjAaj  =  0, 


(4.35) 


where  6  is  the  zero  matrix. 

We  first  contract  Ga  to  K  in  a  coordinate  free  context.  Let  f  be  the 
subalgebra  spanned  by  {Cj  =  2(^2  —  As)}.  A  complementary  subspace  V  is  that 
space  spanned  by  {Cs  —  A2  +  Az^Ci  -  Ai},  so  that  <(tx  is  the  map: 


ojCi  +  02^2  uaCsi*-— ♦AojCi  +  02^3  d"  AosCa. 


(4.36) 
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Calculating, 


\\m4>l'[4>\Ai,2{A2  -Aa)] 

(4.37) 

= 

lim(^^*[AAi,2>l2] 

(4.38) 

= 

lim<^j’(2A)A2 

(4.39) 

= 

linn  ^_x'2A{(i42  —  A3)  +  (.42  +  A3)} 

(4.40) 

= 

l|m  ^^'A(C2  +  C3) 

(4.41) 

= 

Iim2A(C2  +  YC3) 

A— 0  A 

(4.42) 

2C3. 

(4.43) 

Hence:  (Ci,C2]i  =  C3.  All  the  other  structure  constants  can  be  verified  in  the 
same  manner  to  be  0.  Therefore  the  three  dimensional  Lie  algebra  contracts  to  the 
Heisenberg  Lie  algebra. 

We  now  follow  Kalnins  and  Miller  and  repeat  the  above  contraction  in 
coordinates.  We  choose  a  new  path  of  basis  elements  parametized  by  e: 


Ci{t)  =  Cv4i 


/€  0  0\ 
0  0  0 
(0  0  0> 


(4.44) 
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(4.45) 


(4.46) 


[Ci(c),C2(c)]  =  cC2(c)  +  C3(£); 

(C2(£),C3(€)]  =  (C,(£),C2(£)] 

=  0.  (4.47) 

Letting  e  — >  0,  we  obtain  a  new  Lie  Algebra,  whose  basis  has  commutation 
relations: 


(C,(0),C2(0)1  =  C3(0); 

P2(0),C3(0)1  =  (C,(0),C2(0)] 

*  0.  (4.48) 


Again,  we  observe  that  the  affine  Lie  algebra  contracts  to  the  Heisenberg 
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R3  X  R3  — ^  R3 

I  0  X  ^  I  ^ 

X  G\  —*  G\ 

Figure  4.1:  Calculation  of  new  operation  of  the  general  affine  group. 

Lie  algebra.  The  coordinate  approach  will  enable  us  to  work  with  the  affine  Lie 
group  as  the  next  section  will  show. 

4.4.2  Contractions  on  the  Lie  Group  Level 

We  will  now  exponentiate  to  obtain  the  sequence  of  the  corresponding 
Lie  groups  Let  AT  be  a  member  of  Q^.  Since  A'  can  be  expressed  as  A  = 
—PC\  +  oCj  +  0C3,  we  obtain: 


a 

0  \ 

- 1] 

0 

exp( 

0 

1 

0 

)  = 

0 

1 

0 

<  0 

0 

e“-»> 

k  0 

0 

e«c-.  ^ 

This  exponentiation  of  the  Lie  algebra  yields  a  continuum  of  groups  parametrized 
by  e.  The  result  is  a  new  coordinatization  of  G\  associated  to  the  c-Lie  algebra. 
We  would  like  to  know  what  the  multiplication  rule  becomes  in  these  coordinates, 
so  that  we  can  find  the  contracted  rule  when  e— »0.  The  figure  4.1  related  to  3.1 
illustrates  the  process: 


The  map  ^ :  R®  -»  exp((7^)  is  the  map: 


38 


(a,/3,c)^ 


0t 

0  1 


0 


(4.50) 


7 


\  0  0  c 

By  considering  the  Jacobian,  since  <i>  can  also  be  seen  as  a  map  from  to  R^, 
and  applying  L'Hopital’s  rule,  we  can  see  that  ^  is  a  one-to-one  differentiable 
transformation.  Hence  its  inverse  exists: 


:  exp(a:) 

(X  y  0\ 


0  1  0 
VO  0  zj 


R^ 


/yln(^  1,  l.,_,  ,  ,  yln( 


[x-l 


— ln(T),-[ln(«)-|- 


(4.51) 


The  operation,  •,  of  the  group  exp(^‘)  is  simply  the  usual  matrix  multiplication. 
By  seeing  ^  as  a  change  in  coordinates,  we  find  the  corresponding  group  operation 
on  as  follows: 


(a,|3,c)*(a',y9',c')  = 


— fl  -  +  — fe-^'  -  11) 
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il3  + 

C-(/5+/3')  -I  )> 


(4.52) 


Now  let  e  —*  0.  By  expanding  the  exponential  functions  as  power  series,  we 
obtain  the  contracted  multiplication  rule: 


(a,  /3,  c)  ★  (a',  c')  =  (a  +  o',  /3  +  c  +  c'  +  .y— )  (4.53) 

This  group  operation  is  identical  to  the  Heisenberg  group  under  symplectic  coordi¬ 
nates.  Therefore  we  have  contracted  the  affine  group  to  the  Heisenberg  group. 

It  is  instructive  to  consider  the  groups  parametrized  by  c.  We  found  that 
the  group  exp(^*)  is  isomorphic  to  each  exp(^*)  by  an  isomorphism 

(a,/3,c)i — ►(^,c/9,c).  (4.54) 

This  fact  is  easily  verified  by  a  routine  calculation.  We  will  use  this  isomorphism 
in  the  next  chapter. 

4.4.3  Contractions  on  the  Level  of  Group  Representations 

The  group  Ga  has  a  unitary  representation  on  L*{%)  parametrized  by  A  as 
follows  3.3.1: 


(4.55) 
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In  the  previous  sections,  we  always  assumed  implicitly  that  A  =  1,  however,  as  we 
will  see,  the  role  of  A  becomes  more  important  when  contracting  group  representa¬ 
tions. 

To  contract,  we  need  a  representation  A' ^  equivalent  to  above,  obtained 

by  intertwining  by 


6,  (4.56) 

noting  that  any  /  €  L^{Ti)  may  be  expressed  as  g{t)  =  for  some  h{t)  € 

Pin).  Thus: 

A'\a,6,c)<7(<)  =  (4.57) 

o 

Now  we  combine  the  parametrization  <f>  of  in  equation  4.50  with  the  represen¬ 
tation  A'*  with  A  =  A.  Thus  K<  has  the  form: 


<■)»«)  = 
+  ^[1  -  t«*l) 


(4.58) 

(4.59) 


The  proceu  of  letting  e  — »  0  is  so  similar  to  the  other  cases  that  we  will  sim¬ 
ply  state  the  result:  after  contraction,  we  obtain  a  new  representation  A'*’{n,/d,c]: 
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+  a).  (4.60) 

Thus  we  have  obtained  a  representation  of  the  Heisenberg  group,  which  is  what  we 
wanted. 

4.4.4  Contractions  and  Ambiguity  Functions 

We  now  contract  the  ambiguity  functions  themselves.  This  involves  an 
additional  ‘‘trick"  as  will  be  seen  in  this  section. 

As  shown  earlier  3.3.1,  the  wideband  ambiguity  function  is  a  scaled  coeffi¬ 
cient  of  the  group  representation  U^{a,b,c): 

Wj\(c,b)  =  </((),e-'’“'t/"(o,t,c),(()>  (4.61) 

=  0*  r  (4.62) 

J—oo  O 

However,  to  create  our  sequence  of  wideband  ambiguity  functions  to  coverge 
to  the  narrowband  ambiguity  function,  we  use  the  functions  /‘  and  where 
/*  ss  t^f  and  similarly  with  g.  We  also  use  the  <f>  coordinates  for  the  argument 
of  U  and  choose  A  s  and  premultiply  the  ambiguity  function  by  an  appropriate 
factor: 
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e-J«l“fliv,,^.(e-fl‘,±(c-«'-l],„-a)  =  (4.63) 

=  (4.64) 

</,e'¥‘£/t(<-«‘,ile-«‘-ll,€c-a)e‘l“,-  >  (4.65) 


But  the  right  hand  side  of  the  last  bracket  is  just  A  <  as  in  equation  4.56.  So  we 


continue; 


e-W(“f  -  1),  „  -  o)  = 

(4.66) 

<  /,  A-lfe-"',  -  1),«  -  o),  > 

(4.67) 

Letting  c  — »  0,  we  obtain  by  the  previous  section: 


<f,K°(a.fi,c)g>  = 

(4.68) 

<  fit),  e-"'V’"'e'“®y(t  +  o)  >  = 

(4.69) 

f  ^  a)c-*'*^‘dt  = 

V— oo 

(4.70) 

(4.71) 

which  it  essentially  the  narrowband  ambiguity  function,  up  to  phase.  We  have  thus 
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contracted  the  coefficients  of  the  two  groups,  approximating  the  narrowband  am¬ 
biguity  function  by  the  narrowband  via  a  formal  mathematical  apparatus,  without 
any  reference  to  physical  factors. 


Chapter  5 


Volume  and  Contractions 


5.1  Introduction 

We  would  like  to  apply  the  idea  of  contractions  to  gain  some  insight  into 
the  important  concept  of  ambiguity  volume.  After  introducing  ambiguity  volume, 
we  will  explore  the  relationship  between  volume  and  invariant  measures  via  the 
Schur-Godement-  Frobenius  theorem.  We  then  will  study  how  volume  changes  with 
contractions.  Our  aim  is  to  make  rigorous  the  statement  of  Sibul  and  Titlebaum 
(37,  page  86]:  that  the  volume  of  the  wideband  ambiguity  function  is  asymptotically 
conserved  as  the  narrowband  case  is  approached. 

5.2  Ambiguity  Volume 

The  ambiguity  functions  of  interest  in  this  chapter  are  the  auto-ambiguity 
functions,  for  example: 


jv/(*.») = r  + 1)'*’"'*.  (51) 

J—oo 

a  special  case  of  the  cross-ambiguity  function.  These  auto-ambiguity  functions  may 
be  seen  as  functions  of  two  variables  on  their  groups.  The  surface  over  the  x  —  y 
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plane  for  the  narrowband  functions  (x,  y,  Nf{x,y))  called  the  narrowband  ambiguity 
surface. 

Ambiguity  volume,  Vn  is  simply  the  volume  under  the  ambiguity  surface 
squared: 


VM)  =  r  r  \Ni(r,y)\^dxdy.  (5.2) 

Jx^OO  J~oo 

Similarly,  wideband  ambiguity  volume,  Vw{f)  is  defined  to  be: 

Vxy{/)  =  r  r  \W,is,T)f-dr.  (5.3) 

Jo  J—oo  S 

The  interesting  fact  is  that  narrowband  ambiguity  volume  is  always  con¬ 
served: 

Proposition  5.2.1 


Vk(/)  =  r  r  \N,{x,y)\^dxdy  =  ll/lir  (5.4) 

^—00  J-^oo 

that  is,  t/ ll/ll  =  1,  then: 

v«(/)  =  /"  r  \N,(x,y)\^dxdy  =  1.  (5.5) 

J-^oo  J—oo 

Proof  We  reference  the  proof  in  [22,  page  640].  It  also  will  be  a  natural  consequence 
of  the  Godement-Frobenius-Schur  theorem  to  be  stated  later. 
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The  proposition  above  has  been  called  the  ‘‘radar  uncertainty  principle’",  or 
the  “law  of  conservation  of  ambiguity”.  It  can  be  seen  as  saying  that  the  ambiguity 
surface  is  like  a  water  bed:  to  squeeze  a  desired  small  peak  at  a  particular  range 
and  speed  implies  that  other  peaks  will  rise  somewhere.  One  of  the  main  problems 
of  radar  was  to  see  what  kind  of  surfaces  were  ambiguity  surfaces.  One  of  the 
approaches  to  the  solution  of  this  problem  wets  in  seeing  what  kind  of  invariance 
properties  they  exhibited,  which  motivated  Cook  and  Bernfield  [12,  page  68].  The 
radar  conservation  property  was  found  to  be  the  “most  important  ambiguity  func¬ 
tion  constraint  since  it  implies  that  all  signals  are  equally  good  (or  bad)  as  long  as 
they  are  not  compared  against  a  specific  radar  environment”  [12,  page  70].  How¬ 
ever,  the  distribution  of  ambiguity  does  change  from  signal  to  signal.  For  example, 
a  signal  with  short  time  duration  would  produce  a  peak  that  could  resolve  two 
objects  flying  close  together,  but  conservation  would  force  the  objects  velocities  to 
be  unresolvable.  When  coupled  with  the  easily  proven  fact  that  the  maximum  of 
an  autoambiguity  function  occurs  at  the  origin,  the  choice  of  surfaces  becomes  even 
more  limited. 

Wideband  ambiguity  does  not  hold  always,  for: 

Proposition  5.2.2  Let  f  €  Then: 

f  I  |/(u;)|*-.  (5.6) 

Jo  J-^oo  S  Jo  U! 


Proof  We  again  reference  [22,  page  643]  for  details. 
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There  are  certainly  /  €  L^{11)  such  that  /o°°  |/(u>)P^  is  not  defined.  The  loss  of 
conservation  in  the  wideband  realm  causes  difficulty  for  selection  of  surfaces.  Only 
one  constraint  in  wideband  signal  design  remains.  An  approach  to  classifying  all 
wideband  ambiguity  surfaces  is  detailed  in  [1]. 

5.3  Volume  and  Group  Theory 

To  apply  group  theory  for  insight  into  the  volume  problem,  we  need  a  few 
definitions. 

Definition  5.3.1  Let  U  be  a  unitary  representation  of  a  Lie  group  G  on  a  Hilbert 
space  H .  Then  U  is  said  to  be  a  square  integrable  representation  if  there  is  a  vector 
g  €  H  such  that: 


[  I  <  U{x)g,g  >  l’d^(x)  <  oo.  (5.7) 

JG 

Such  a  vector  g  is  called  admissible. 

We  now  quote  without  proof  the  Froebenius-Schur-Godement  theorem: 

Proposition  5.3.1  If  U  is  a  square  integrable  representation  of  a  group  G  on  a 
Hilbert  space  H,  there  exists  a  unique  self-adjoint  operator  Q  such  that: 

i)  The  set  of  admissible  vectors  coincide  xoith  the  domain  of  Q. 

ii)  Let  gi  and  g2  be  two  admissible  vectors  and  let  /]  and  /j  be  any  vectors 


in  H.  Then: 


48 


<  f\,U{x)gi  >  <  f2,U{x)g2  >dn(x)  =<  Qg2,Qg\  ><  /i,/?  >  •  (5.8) 

Hi)  Finally,  ifGis  unimodular,  then  Q  is  a  multiple  of  the  identity  operator; 
hence  the  domain  of  Q  will  be  all  of  H,  that  is,  all  of  H  constitutes  the  set  of 
admissible  vectors  for  the  representation. 

The  proposition  is  proved  in  detail  in  (38,  page  62]. 

Now  we  will  apply  these  new  facts  to  our  ambiguity  functions.  The  equa¬ 
tion  5.7  may  be  recast  as  stating  that  U  is  square-integrable  if  there  is  a  signal  g 
such  that  its  ambiguity  volume  is  finite.  This  general  statement  applies  to  both 
narrow  and  wideband  when  considering  the  appropriate  representation  as  indicated 
previously.  The  equations  5.5  and  5.6  show  that  under  this  formulation,  (/«  and 
Ua2  arc  square  integrable  representations  of  their  respective  groups. 

Since  the  Heisenberg  group  is  unimodular  [22,  page  640],  all  of  L^iTl)  is 
admissible  by  the  Frobenius-Schur-Godement  theorem.  This  is  not  true  for  the 
affine  group:  there  are  functions  for  which  the  equation  5.6  will  not  be  finite  on  the 
right  hand  side.  And  the  Frobenius-Schur-Godement  theorem  does  not  guarantee 
that  all  of  would  be  admissible  for  .  By  the  converse  of  the  third  part 

of  the  theorem,  we  see  that  the  affine  group  cannot  be  unimodular.  And  sure 
enough,  it  is  not  [22,  page  643].  We  can  see  that  there  is  some  connection  between 
unimodularity  of  groups  and  ambiguity  volume. 
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We  wanted  to  see  how  invariant  measures  behaved  in  the  e-groups  involved 
in  the  contraction  of  the  affine  group  to  the  Heisenberg  group. 

Proposition  5.3.2  The  left  invariant  measure  of  an  e-group  is 

1  - 

— - - dad0dc  (5.9) 

pe 

and  the  right: 


-  1 

— - - dad^dc. 


(5.10) 


Proof  We  may  prove  the  proposition  by  straightforward  verification.  For  example, 
the  right  invariant  measure  may  be  verified  by  showing  that: 


j  /((a,  /S,  c)  •  (o',  6',  c')]dadl3dc  = 

(5.11) 

f  f{u,v,w)^  dudvdw  = 

J  o{u^v,w) 

(5.12) 

J  f{u,v,w)dudvdw, 

(5.13) 

where  u,  v,  and  w  denote  the  components  of  the  right  hand  side  of  equation  4.52 
in  their  respective  orders,  and  denotes  the  Jacobian  of  the  transformation, 

computed  by  taking  the  reciprocal  of 
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An  interesting  aside  is  that  the  isomorphism  from  the  previous  chapter, 
carries  the  left  invariant  measure  of  the  group  when  c  =  1  to  the  left  invariant 
measure  of  any  specified  c-group. 

As  long  as  c  >  0,  we  can  see  that  the  groups  will  never  be  unimodular, 
but  as  c  — »  0,  the  groups  become  more  and  more  unimodular.  It  seems  that  the 
Frobenius-Schur-Godement  theorem  might  help  us  to  see  that  ambiguity  volume 
becomes  conserved. 

Thus  the  statement  that  the  “wideband  ambiguity  volume  is  asymptotically 
conserved”  now  has  a  mathematical  foundation. 


Chapter  6 


Conclusions  and  Questions  for  Further  Research 

We  have  shown  how  group  theory  can  provide  insight  into  some  of  some 
problems  in  RADAR  and  SONAR  signal  processing.  We  began  by  noting  that 
the  RADAR  and  SONAR  ambiguity  functions  are  coefficients  of  the  unitary  rep¬ 
resentations  of  the  Heisenberg  and  affine  groups  respectively.  This  basic  fact  not 
only  unifies  the  various  forms  of  the  wideband  ambiguity  function,  but  helps  us  in 
simplifying  various  computations.  We  also  showed  how  contractions  helped  model 
the  approximations  of  the  wideband  to  the  narrowband  ambiguity  function  under 
suitable  conditions.  And  we  showed  that  the  wideband  ambiguity  volume  is  not 
conserved,  even  when  the  affine  group  in  question  is  very  close  to  the  Heisenberg 
group.  Finally,  we  showed  that  harmonic  analysis  on  finite  groups  may  provide 
insight  into  creating  another  transform. 

Our  work  is  just  the  beginning  of  a  '^wideband”  version  of  Schempp's  work 
[38].  For  example,  the  characterization  of  narrowband  ambiguity  surfaces  was 
shown  to  be  a  natural  consequence  of  harmonic  analysis  on  the  Heisenberg  group 
Schempp.  He  also  solved  the  invariant  problem,  finding  the  group  under  whose 
motions  the  narrowband  ambiguity  surfaces  remuned  invariant.  A  wideband  coun¬ 
terpart  to  these  results  seems  daunting,  since  the  affine  group  is  not  nilpotent  as 
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the  Heisenberg,  but  only  solvable.  To  complicate  Schempp’s  approach  further,  the 
af&ne  group  is  not  unimodular.  The  theory  of  positive  definite  functions  on  the 
affine  group,  or  the  autoambiguity  functions  in  this  thesis,  is  not  very  well  known 
for  that  reason. 

We  also  would  have  liked  to  characterize  the  range  of  the  wideband  autoam¬ 
biguity  function  as  Auslander  and  Tolimieri  did  with  the  narrowband  autoambigu¬ 
ity  function  [4].  However,  the  affine  group  seems  immune  to  their  approach,  since 
it  lacks  a  discrete  subgroup  to  quotient  it  to  a  compact  manifold,  to  enable  us  to 
apply  some  theory  of  “solv-theta”  functions  to  it  as  they  did. 

Many  questions  remain  to  be  answered.  For  example,  what  is  the  geometric 
picture  of  this  contraction?  Dooley  shows  some  elegant  geometric  pictures  for 
some  of  semisimple  group  contractions  [15].  Can  Kirillov  coadjoint  orbital  theory 
be  used  to  show  what  is  happening  with  our  two  groups?  Here,  the  difficulty 
is  that  the  groups  concerned  lack  “semisimplicity”,  and  are  solvable.  We  have 
calculated  the  coadjoint  action  of  the  e-groups,  and  the  action  seems  to  *"biow 
up'’  at  the  limit.  What  this  means  we  are  not  sure.  Is  the  relationship  between 
wideband  and  narrowband  analogous  to  the  relationship  between  quantum  and 
classical  mechanics?  Is  there  a  method  of  ‘‘quantization”  from  a  signal  processing 
point  of  view? 

The  discrete  affine  group  harmonic  analysis  also  provides  many  interesting 
questions.  Is  there  a  finite  cross  ambiguity  function?  How  should  the  signals  be 
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sampled?  Is  there  a  solvable  analog  to  Shannon’s  sampling  theorem?  What  is  the 
significance  of  the  afRne  Fourier  transform?  The  affine  group  here  seems  to  lack  the 
beautiful  “Stone-von  Newman”  property  of  the  Heisenberg:  that  the  Heisenberg 
has  essentially  one  unitary  representation. 

Finally,  we  hope  to  see  more  areas  of  signal  processing  “invaded”  by  modern 
mathematics.  While  the  insights  gained  would  be  invaluable,  perhaps  new  devices 
to  serve  the  needs  of  humanity  w’ill  be  created  as  a  result. 


Appendix 


The  Discrete  Affine  Fourier  Transform 

A.l  E^TRODUCTION 

Discrete  Fourier  Transforms  (DFT)  have  been  standard  tools  in  digital 
signal  processing  for  approximately  three  decades.  If  we  view  the  theory  of  the 
Discrete  Fourier  Transform  as  harmonic  analysis  on  a  group  of  integers  modulo  a 
prime  number,  we  can  gain  new  insights  that  will  allow  us  to  develop  new  discrete 
Fourier  transforms  as  harmonic  analysis  on  other  discrete  groups. 

Group  theory  enables  us  to  see  various  transforms  in  a  unifying  manner. 
For  example,  the  Fourier  transform  is  intimately  related  to  the  group  of  the  real 
numbers;  the  Z-transform  is  connected  with  the  integers;  and  the  Discrete  Fourier 
Transform  is  related  to  the  group  of  integers  mod  n.  We  introduce  a  new  transform 
by  considering  the  finite  affine  group,  and  find  that  all  the  theory  of  the  Discrete 
Fourier  Transforms  generalizes  in  a  pleasant  way.  The  formula  for  the  Discrete 
Affine  Transform  remains  essentially  the  same  as  the  Discrete  Fourier  Transform. 

The  underlying  theory  for  this  transform, namely  that  of  harmonic  analysis 
on  finite  groups,  has  been  known  for  quite  some  time  [23,  24].  We  have  essen¬ 
tially  applied  the  theory  to  the  finite  affine  group,  having  been  motivated  by  radar 
ambiguity  function  theory  [36,  9,  5). 
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We  will  review  familiar  properties  of  DFT’s  that  are  consequences  of  har¬ 
monic  analysis  on  finite  fields,  auid  then  investigate  harmonic  analysis  on  the  finite 
affine  group  Ap,  where  p  is  a  prime  number.  Focusing  on  a  specific  finite  affine 
group,  we  then  develop  the  Discrete  Affine  Fourier  Transform  (DAFT),  using  a 
specific  finite  affine  group  for  illustration.  (This  transform  should  not  be  confused 
with  the  Affine  Wavelet  Transform).  After  establishing  the  Inverse  Discrete  Affine 
Fourier  Transform  (IDAFT),  some  transform  properties  analogous  to  the  DFT  will 
be  shown,  namely  those  involving  the  relationships  between  the  DAFT  and  the 
norms  and  convolutions  of  its  function  spaces.  Finally,  we  show  how  the  DAFT 
would  be  useful  for  solving  wideband  inverse  problems. 

A.2  Some  More  Mathematical  Background 

Associated  with  a  set  5  is  the  vector  space  of  all  complex  valued  functions 
on  5  denoted  by  ^{S).  This  T{S)  is  a  vector  space  in  a  natural  way.  If  5  is  finite, 
^(S)  is  isomorphic  to  where  n  is  the  number  of  objects  in  5. 

As  an  example,  ^{Zn)  is  the  space  of  all  complex  valued  sequences  on 
=  (/o»*  • -i/n-i)}-  Addition  on  ^(Zn)  is  simply  componentwise  addition 
of  sequences  and  scalar  multiplication  is  defined  similarly.  An  norm  may  be  defined 
on  ^(Zn):  ll/ll*  =  /*?*'  These  norms  induce  an  inner  product  as  in  usual 

Hilbert  space  theory. 

A  * 

The  space  ^{G)  is  defined  analogously,  with  G  forming  the  index  set  of  the 
sequences.  Each  member  pk  of  the  sequence  will  be  a  matrix  of  fixed  dimension 


rik.  A  norm  for  /  =  is:  |!/1P  =  E)r=o”t<r(//*),  where  ’  denotes 

the  Hermitian  transpose  of  a  matrix.  In  the  case  of  .2„,  when  all  the  irreducible 
representations  are  one  dimensional,  the  norm  becomes  the  familiar  one.  Of  course, 
these  norms  induce  their  corresponding  inner  products. 

A.3  The  Discrete  Fourier  Transform 

As  stated  above,  the  theorj'  of  the  Discrete  Fourier  Transform  (DFT)  is 
merely  harmonic  analysis  on  Zn^  the  group  of  integers  modulo  an  integer  [28,  pages 
546-548].  Harmonic  analysis  simply  means  expressing  the  functions  on  the  group 
as  a  sum  of  the  irreducible  representations  of  the  group. 

The  DFT  can  be  seen  as  the  mapping  from  the  space  of  all  complex  se¬ 
quences  on  Z„  to  the  space  of  complex  sequences  on  Z„.  Its  explicit  form  is: 

(A.i) 

(sO 

The  DFT  formula  is  thus  intimately  related  to  the  ‘^harmonics”  of  Z„. 

In  addition,  the  Inverse  Discrete  Fourier  Transform  (IDFT)  may  be  defined,  also 
involving  the  characters  of  Zn- 

=  (A'2) 

P  wsO 

The  DFT  has  some  nice  properties.  It  is  well  known  that  the  DFT  preserves 
norms:  |(/|I^  =  ^||/||*-  Secondly,  the  (cyclic)  convolution  of  two  n-length  sequences 
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defined  by: 


if*9)t  =  '^fT9t-r,  (A.3) 

T=1 

t  —  T  performed  with  modulo  n  arithmetic,  is  conveniently  mapped  to  component¬ 
wise  multiplication  of  their  transforms:  {f*g)^  =  /w  •  These  two  properties  of 
the  DFT  form  a  foundation  for  a  full  range  of  applications  [33]. 

A.4  The  Discrete  Affine  Fourier  Transform 

A.4.1  The  Finite  Affine  Group  and  Its  Representation  Theory 

Motivated  by  the  affine  group’s  intimate  relationship  to  the  wideband  am¬ 
biguity  function  [36,  9],  and  inspired  by  Auslander  and  Tolimieri’s  approach  to  the 
DFT  [5,  page  866],  we  investigate  the  so  called  finite  affine  group,  denoted  by  Ap. 
This  is  the  set  of  2  by  2  matrices  of  the  form: 


s 

0 


(A.4) 


with  s  taking  nonzero  values  in  Zp  and  r,  any  values  in  Zp.  The  group  operation  will 
be  the  usual  matrix  multiplication,  however  using  modulo  p  arithmetic.  Modulo  p 
arithmetic  is  necessary  for  every  nonzero  integer  to  have  a  multiplicative  inverse. 
This  condition  may  be  loosened  slightly  by  requiring  that  the  multiplicative  group 
of  the  ring  have  a  primitive  generator,  a  condition  true  for  £„  with  n  =  2, 4,  p“ ,  2p° 
[25,  page  44).  However,  for  simplicity’s  sake,  we  shall  consider  only  the  case  of  Ap. 
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A. 4. 2  The  Discrete  Affine  TVansform 

The  theory  of  Fourier  transforms  for  the  finite  noncommutative  group  had 
been  discovered  for  quite  some  time  (24,  page  114]  .  We  will  now  construct  the 
Discrete  Affine  Wavelet  Transform  (DAFT),  the  affine  analogue  of  the  DFT,  using 
this  theory.  This  process  involves  identifying  “harmonics”  on  the  group  Ap,  and 
then  applying  the  formula  for  the  Fourier  transform  in  the  literature  to  this  case 
[24,  pages  77,81).  We  shall  use  A3  as  an  illustrative  example. 

Identifying  all  the  harmonics  of  Ap. 

We  identified  all  the  irreducible  representations  of  Ap  in  a  semi  ad  hoc 
manner,  by  using  two  results  in  finite  group  representation  theory.  The  first  is  that 
the  number  of  the  characters  of  a  group  is  equal  to  the  number  of  its  conjugacy 
classes,  and  the  second  is  that  the  sum  of  the  squares  of  the  dimensions  of  each 
irreducible  representation  is  equal  to  the  number  of  elements  of  the  group  [29].  We 
found  that  there  were  p  conjugacy  classes  in  Ap,  hence  p  irreducible  representations. 
And  then: 


1^4-...  +  1^.4-(p-  1)^  =p(p~  1).  (A.5) 

p-ittmes 

Therefore  Ap  had  to  have  p  one  dimensional  characters  and  one  p  —  1- 
dimensional  character.  We  found  the  explicit  form  of  the  representations  by  some 
trial  and  error,  by  consideration  of  their  subgroups  with  some  trial  and  error. 
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Incidentally,  we  found  [34],  which  confirmed  some  of  the  results  of  this  section. 
The  explicit  form  of  the  representations  follow.  VVe  note  that  not  all  of  them  are 
one  dimensional. 

The  "harmonics”  of  Ap  consists  of  p  —  1  one  dimensional  representations 
and  one  p  —  l*dimensional  representation,  a  total  of  p  elements,  as  follows: 


Po(fif',T)  =  1 

Pt(g  ,r)  =  cp-> 


t) 


2ntkM 

e 


Pp-7(ff\r)  =  c  »>-»  (A.6) 

where  p  is  the  number  with  which  every  nonzero  element  of  Zp  may  be  expressed 
as  a  power  of  g.  The  p  —  1-dimensional  representation  Pp.i  is  that  representation 
induced  from  the  one  dimensional  representation: 

x(l,T)  =  e*V^  (A.7) 


The  difficulty  of  expressing  the  DAFT  explicitly  in  general  results  from  not  yet 
knowing  the  explicit  form  of  this  representation. 
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We  now  consider  the  representation  theory  of  >^3.  The  group  A3  consists 
of  six  elements,  since  there  are  two  choices  for  s  and  three  choices  for  t.  Because 
s  is  never  zero,  we  can  express  it  as  a  power  of  2.  (This  fact  holds  for  general  Ap, 
although  the  base  of  the  exponential  may  not  necessarily  be  2).  Hence  the  space 
of  complex  functions  on  the  group  are  those  sequences  of  length  six.  However, 
the  set  of  “harmonics”,  the  irreducible  complex  representations  of  .43,  consists  of 
only  three  elements,  a  consequence  of  the  group’s  noncommutativity.  Hence  we  no 
longer  have  a  group  of  harmonics,  but  something  more  complicated. 

Explicitly,  the  characters  of  A3  are: 


(A.8) 


Pi(2%r) 


(A.9) 

(A.IO) 


Constructing  the  DAFT 

The  next  two  sections  consist  of  applying  the  results  of  Hewitt  and  Ross  [23, 
24]  to  our  specific  group,  involving  translation  of  their  mathematical  formulations. 
Given  a  sequence  /  on  Ap,  that  is,  /  has  the  form: 


/  —  (fifif  •  •  • » fi,pj  /3.O1  •  •  • » fp-i,p)t 


(A.ll) 


*  «k 

we  can  construct  a  sequence  /  on  Ap,  the  set  of  “harmonics"  of  Ap,  called 


the  Discrete  Fourier  Transform  (DAFT)  as  follows: 

Since  each  representation  of  Ap  is  associated  to  a  finite  linear  space,  we 
can  use  the  space’s  inner  product  to  define  a  coefficient  of  two  vectors,  u  and  v, 
associated  with  the  representation: 

C^y{9\T)  =<  /)(5*,t)u,v  >  .  (A. 12) 

We  use  this  concept  to  clarify  the  transform  for  a  finite,  noncommutative  group 
such  as  ours.  The  DAFT  is  simply  the  sum: 

<  /(«)U,V  >=  Y.  Q,/(s‘,r)  (A.13) 

* 

For  Ap,  fiPk)  can  be  seen  as  an  n  x  n  matrix,  where  n  is  the  dimension  of  the  linear 
space  involved. 

The  DAFT  may  be  written  explicitly  for  A3  as  follows:  if  /  =  (/lo,  /n,  /12,  /20,  /211  /22) 
is  a  sequence  on  .As,  the  sequence  /  consists  of: 

fpo  =  JZ  /2*.t 

(2*.t)€^^ 

L  =  E  (-l)V2-.r 

(2*,t)€-4, 

* 

/«  = 

/lO  +  +  W^/12  /20  +  ‘*7^/21  +  ‘*'/22 

/20  +  w/ai  +  u;*/22  fio  +  w^/ii  +  w/ij 
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Constructing  the  IDAFT 

The  DAFT  is  not  onto  for  general  groups.  However,  since  Aj,  is  finite,  the 
DAFT  is  onto,  and  we  have  the  existence  of  the  Inverse  Discrete  Affine  Fourier 
Transform  [24,  page  420].  Consider  the  space  of  all  matrix  valued  sequences  on  Ap, 
that  is,  the  space  of  all  sequences  with  the  “harmonics”  of  Ap  as  the  index  set: 


(A.14) 


where  the  first  p  ~  1  entries  are  simply  scalars,  and  the  last  entry  isap— 1  xp— 1 
matrix.  The  function  /  whose  DAFT  is  /  is: 


ib=0 


+(p~  l)tr(/fc-pp.j(p*,r)),  (A.15) 


where  tr  denotes  the  trace  of  a  matrix.  This  formula  will  be  called  the  Inverse 
Discrete  Affine  Fourier  Transform  (IDAFT). 

The  IDAFT  can  be  written  for  the  case  of  A3  as  follows;  Let  /  =  {/o,  /i,  A) 
be  a  sequence  on  the  set  of  “harmonics”  of  A3.  We  note  in  passing  that  /a  is  a  two 
by  two  complex  matrix.  Then  the  function  whose  DAFT  is  /  is: 


/j'.T  =  /o  +  (—I)*/!  +  2tr(/2  •  P2(2*,t)) 


(A.16) 
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A. 4.3  Interesting  Properties  of  the  DAFT 

The  DAFT  has  a  number  of  interesting  properties,  quite  similar  to  the 
DFT.  We  will  mention  two  here,  one  related  to  norms  and  the  other  related  to  con¬ 
volutions.  We  have  proved  them  explicitly  for  .43  only,  however  the  generalizations 
seem  straight  forward.  Since  the  proofs  are  computational  in  nature,  they  have 
been  omitted. 

Defining  a  norm  on  the  sequences  on  A3  by; 

ll/ll  ==  l/ol  +I/1I  +2fr(/3/3),  where  *  denotes  the  Hermitian  transpose,  we 
found  that:  ||/||  =  61I/||^.  This  fact  should  generalize  to  be:  ||/||  =  p{p  —  l)a'/||^, 
for  general  Ap. 

The  afiine  group  theoretical  convolution  of  the  form: 

f*g{s,T)=  Y,  fir^y)9[{x,y)~Hs,T)]  (A.17) 

(*,»)€.4p 

is  mapped  via  the  DAFT  to  componentwise  multiplication: 

Fg  =  H.  (A.18) 

The  right  side  denotes  componentwise  multiplication,  as  in  the  DFT  case, 
with  the  matrix  components  multiplied  as  regular  matrices.  We  were  surprised  that 
the  products  reversed  order  under  the  transform,  a  fact  proven  for  the  case  of  A3 
only  (This  fact  has  not  been  mentioned  in  any  works  that  we  have  seen). 
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A.5  Towards  An  Application  of  the  DAFT 

This  section  will  describe  a  possible  application  of  the  ideas  discussed.  It 
will  be  seen  that  it  can  only  be  a  sketch  at  best,  yet  hopefully  detailed  enough  to 
motivate  further  research. 

One  of  the  important  applications  of  affine  deconvolution  is  the  solution  of 
wideband  inverse  problems.  Consider  the  received  signal  r(<)  given  by; 

r(l)  =  /  /  c"“p(so,To)/(e“'“t  -  To)dTodso  (A. 19) 

where  p{sq,tqi  is  the  complex  scattering  strength  at  delay  Tq  and  dilation  so,  and 
/(<)  is  the  transmitted  probing  signal.  The  time  varying  transmission  or  scattering 
channel  causes  the  signal  to  be  delayed  by  Tq  and  its  time-scale  dilated  by  e"^. 
For  narrowband  signals  the  time  dilation  reduces  to  the  well  known  Doppler  shift. 
The  inverse  problem  is  to  determine  p(5o, ^o)  from  received  signal  r{t)  and  from 
the  known  probing  signal  f{i).  This  can  also  be  called  the  identification  of  linear 
time- varying  systems  [46]. 

Let  us  consider  processing  of  the  received  signal  by  the  correlation  receiver. 
The  correlation  receiver  computes  an  inner  product  of  received  signal  r(t)  and  the 
processing  waveform  g{i),  which  is  time-delayed  and  scaled.  If  the  received  signal 
contuns  colored  noise  in  addition  to  the  scattered  probing  signal,  the  optimum 
processing  waveform  is  not  the  dilated  and  delayed  version  of  the  probing  signal. 
We  also  recognize  that  the  correlation  receiver  output  /(s,r)  is  an  affine  wavelet 
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transform  of  r{t)  with  respect  to  the  processing  waveform  g(t).  Thus, 


/(s,t)  =  Wgr  (A.20) 

=  <r{t),g(e~*t-T)e~"  >  (A.21) 

g’{e~*t  —  T)e~^dTQdsQ^  dt  (A. 23) 


Interchanging  the  order  of  integration,  we  have: 

f  f  (I  -  To)  g'’{e~i)e~‘dt)dTodso.  (A.24) 

J—ooJ  0  \J—oo 

We  now  recognize  that  the  inner  integrals  are  an  inner  product  and  we  have: 

/  pisoyTo)  <  U{so,To)f{t),U{s,T)g{i)  >  drodso,  (A.25) 
J—ooJ  0 

where  U{so,to)  and  U{s,t)  are  unitary  affine  transformations  [9].  Since  U{so,to) 
is  unitary,  its  adjoint  is  equal  to  its  inverse  and  we  have: 


/(s,r) 


/  /  p(sq,To)  <  f{i),[/'^{so,To)Uis,T)g(t)  >  dTodao  (A.26) 
J—e»J  0 

/  /  p(«o,Tb)  {A.27) 

J—ooJ  0 
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(/~  f(t)  +  c'^*"*** Vo  +  T)e-^dt  )dTodso  (A.28) 

The  inner  integral  is  the  wideband  ambiguity  function, 

W Hence  the  most  recent  equation  may  be  seen  as  a 
continuous  convolution  on  the  affine  group  [32,  6]: 

l{s,T)  =  p*WAf,g  (A.29) 

=  /  /  p{sq,to)W Ajj,[{t\T){e‘°,To)~']dTodso  (A.30) 

J—ooJ  0 

=  /  /  -  c~‘*“'®Vo)dTodso)-  (A.31) 

Thus  the  solution  of  the  wideband  inverse  problem  is  equivalent  to  deconvolution. 

The  solution  of  this  inverse  problem  can  be  attacked  in  several  ways.  One 
could  use  Kahlil’s  formula  as  suggested  by  Napharst  [32],  or  wavelet  transform 
inversion  as  suggested  by  Heil  and  Walnut  [22],  Young  [46],  and  Grossmann  and 
Morlet  [17].  Use  of  the  DAFT  to  deconvolve  a  discrete  scattering  problem  would 
be  another  approach  for  additional  insight  into  the  solution  of  the  inverse  problem. 
Many  difficult  computational  and  discretization  problems  remain  to  be  worked 
out.  We  hope  that  the  Discrete  Affine  Fourier  Transform  Theory  may  have  other 
interesting  applications  to  wideband  time-varying  system  theory. 
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A.6  Looking  Ahead 

We  have  performed  some  harmonic  analysis  on  a  finite  noncommutative 
group.  Noncommutativity  generally  makes  such  analysis  very  difficult,  but  finite¬ 
ness  of  the  groups  saves  the  day.  The  same  strategy  may  be  applied  to  other 
groups  such  as  the  Heisenberg  group  for  instance.  As  the  affine  group  underlies 
the  wideband  ambiguity  radar  theory,  the  Heisenberg  group  plays  a  major  role  in 
narrowband  ambiguity  radar  theory.  Deconvolution  is  a  major  problem  there  too. 
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